Universal dephasing rate due to diluted Kondo impurities 
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We calculate the dephasing rate due to magnetic impurities in a weakly disordered metal as 
measured in a weak localization experiment. If the density ns of magnetic impurities is sufficiently 
low, the dephasing rate l/r v is a universal function, \/t v = (ns / 'v) f(T /Tk) , where Tk is the Kondo 
temperature and v the density of states. We show that inelastic vertex corrections with a typical 
energy transfer AE are suppressed by powers of 1/(t v AE) oc ns. Therefore the dephasing rate 
can be calculated from the inelastic cross section proportional to TrulraT — 7tj/T| 2 , where T is the 
T-matrix which is evaluated numerically exactly using the numerical renorrnalization group. 



Dephasing, i.e. the loss of wave coherence, is a ubiqui- 
tous phenomenon in the quantum mechanics of complex 
systems. It is of relevance to any experiment where both 
interference and interactions play a role and is, there- 
fore, of profound importance in all areas of nanoscopic 
and mesoscopic physics. 

While the basic phenomenon of dephasing as such is of 
rather general nature, the concrete definition of a dephas- 
ing rate, and its experimental determination vary from 
context to context. In this Letter we consider the dephas- 
ing rate as determined by weak-localization (WL) mea- 
surements in metals Q. In weakly disordered metals, the 
interference of electronic wave functions on time-reversed 
paths leads to a characteristic reduction (or enhancement 
in the presence of spin-orbit interactions) of the conduc- 
tivity. The magnitude of this effect is controlled by the 
dephasing time t v — the typical time-scale over which 
electrons get entangled with their environment (phonons, 
other electrons or dynamical impurities) thereby loosing 
the ability to interfere. Even small magnetic fields break 
time-reversal invariance thus prohibiting the interference 
of time-reversed paths. Fitting the magnetoresistivity to 
WL theory is a means to determine the dephasing rate 
1/t v with high precision. 

Surprisingly, most of these experiments [2| show a sat- 
uration of the dephasing rate at the lowest accessible 
temperatures T, while theoretically it is expected that 
in the limit T — > all inelastic processes freeze out when 
the system approaches its (time-reversal invariant and 
unique) ground-state. This has lead to an intense dis- 
cussion [3, 0, 0, HI as to whether quantum fluctuations 
can induce dephasing at T = 0. While we believe that 
this latter scenario is theoretically excluded for electrons 
in a disordered metal, the presence of only a few parts- 
per-million (ppm) of dynamical impurities — realized by 
atomic two-level systems or by magnetic impurities 
IE EJ EH — may be an alternative cause of the satura- 
tion phenomenon. Indeed it has been shown experimen- 
tally that magnetic impurities lead to an apparent satu- 



seem to follow the predictions of Altshuler, Aronov and 
Khmelnitsky for the dephasing induced by Coulomb 
interactions. 

The effect of dynamic magnetic impurities on 1 /t v has 
first been considered by Ohkawa, Fukuyama and Yosida 
|l2| (for the static case see Ref. 0) using perturbation 
theory (generalized to renormalized perturbation theory 
in Refs. 0, Q|) which limits the range of applicability 
to temperatures T 3> Tk larger than the Kondo tem- 
perature, Tk- For T <C Tk, a quadratic T dependence, 
1/t„ oc T 2 has been predicted based on Fermi liq- 
uid arguments. In this paper we argue that the impact 
of static and diluted dynamic impurities can largely be 
treated separately. This separation enables us to com- 
bine a perturbative approach to the static disorder with 
a (numerically) exact treatment of the Kondo interac- 
tion. In this way, the effect of a small concentration, ns, 
of magnetic impurities on 1/t v can be explored in the 
entire crossover range from T 3> Tk down to T <C Tk- 

We consider the Hamiltonian H = Hq + Hg , where 



d rf x4(x) |--n + V(x) c CT (x) (1) 



describes electrons in a weak non-magnetic disorder po- 
tential V, modeled for convenience by Gaussian white 
noise: (V(x)V(x'))v — 6(x — x') where v is the den- 
sity of states (DOS), r = l/vp the elastic scattering time 
and I the elastic mean free path. (None of our results is 
affected by the precise choice of band-structure or model 
of disorder.) The coupling of the electrons to a small 
concentration ng of spin-1/2 impurities is described by 
the Kondo Hamiltonian 



H s = J J~] Sj{4.(xi)cw cy (xj 



(2) 



To calculate physical quantities, we have (i) to average 
over V(x) and (ii) the positions x, of the spins, taking 
into account (iii) the exchange coupling J to all orders. 
Two small parameters will help us in attacking this prob- 



lem: The concentration of magnetic impurities ns is tiny 
ration of 1/t v at least in some T range [9|. In contrast, [more precisely us/tyTa) <C 1] and also the ratio of elec- 
tronic wave-length and elastic mean free path, l/(kpl), 
is small but finite. 



some extremely pure Au and Ag samples with negligible 
concentration of impurities |lCj show no saturation and 
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FIG. 1: Bethe-Salpeter equation for the Cooperon C in the 
presence of (dilute) magnetic impurities to linear order in ns . 
Co is the bare Cooperon in the absence of interactions and V 
the irreducible vertex obtained by adding self-energy, elastic- 
and inelastic vertex contributions. The crosses with attached 
dashed lines denote the averaging over impurity positions X;, 
the squares the interactions to arbitrary order in J. 



We wish to explore the impact of dynamic impu- 
rity scattering on the WL corrections to the electric 
conductivity. Technically, this amounts to computing 
the impurity generated 'self energy' or 'mass' of the 
Cooperon describing the coherent propagation of an elec- 
tron and a time-reversed hole in the disordered environ- 
ment. As justified in detail below, we neglect mixed inter- 
action/disorder diagrams. The problem therefore reduces 
to the solution of the Bethe-Salpeter equation depicted 
diagrammatically in Fig.^a). To linear order in ns, the 
two-particle irreducible vertex T shown in Fig. ^b) can 
be separated into three distinct contributions: self-energy 
diagrams (the first two terms in Fig. ^b)), an 'elastic' 
vertex correction with no energy transfer between upper 
and lower line, and an 'inelastic' vertex where interac- 
tion lines connect the two lines. We begin by focusing on 
the elastic contributions to the Cooperon self energy as 
inelastic contributions give vanishingly small corrections 
for small ns, see below. 

Since self-energy and elastic vertex contributions con- 
serve the energy of single electron lines, the solution of 
the reduced Bethe-Salpeter amounts to a straightforward 
summation of a geometric series. Setting the center-of- 
mass frequency SI (see Fig.QJ to 0, the Cooperon obtains 
as C« =0 (e,e>) = ^S(e - e')/{Dcf + l/r v {e,T)) and 
the WL correction to the conductivity is given by 
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with the T and e dependent dephasing rate 

1 = ^[nvIm[T A (e)]~\nvT R (e)\ 2 ] 



T v {e,T) ttv 



(3) 



(4) 



The many-body T-matrix [defined by the Green func- 
tion G xx <(e) - G xx ,(e) + G x0 (e)T(e)G° x , (e)] describes 
the scattering of electrons from a single magnetic impu- 
rity. Eq. 0} affords a simple interpretation: Im[T A (e)] 
is proportional to the total cross section, while | T R (e)| 2 
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FIG. 2: a) Interference of interacting electrons, b) Cor- 
rections due to the inelastic vertex to lowest order, c) WL 
correction to the dephasing rate linear in ns- d) Altshuler- 
Aronov-Khmelnitsky type corrections arising to order n§. 



describes the elastic cross section. Their difference is - 
by definition - proportional to the inelastic cross section 
recently been introduced by Zarand et al. |l5j . Note that 
the optical theorem guarantees the vanishing of the in- 
elastic cross-section for static impurities. Together with 
the numerical evaluation of 1/t v and the analytic esti- 
mates of the leading corrections, given below, Eq. is 
the main result of this paper. 

To determine the range of applicability of Eq. 0}, 
three classes of corrections have to be considered: con- 
tributions from the inelastic vertex, mixed interaction- 
disorder diagrams, and higher-order corrections in rig- 
We begin by considering the lowest order expansion of 
the WL correction in the inelastic vertex (cf. Fig.^b)), 



de / dSl f(e, fi)r in (e, Sl)I d (e, Si). (5) 



Here, / denotes some thermal function restricting e and SI 
to values smaller than T, Lj n is the inelastic vertex and 

Id(e,Si) ~ [(l/r^) 2 + fi 2 ] 4 results from momentum- 
integration over two Cooperons in d dimensions. For suf- 
ficiently small ns <C vTk the typical energy AE ~ SI ~ T 
[lfij | exchanged at the inelastic vertex greatly exceeds the 
dephasing rate (see below) 1/t v oc ns- Evaluating the 
integral © for AE 3> we obtain 



Aai ( 1 



AdWL V AE Tu; 



(4-d)/2 



< 



(4-d)/2 



«1. (6) 



It is straightforward to verify that this estimate pertains 
to higher orders in the expansion in the inelastic vertex. 

The irrelevancy of the inelastic vertex for processes 
with typical energy transfer larger than 1 /t m has been ob- 
served before by several authors (e.g. [ll|,llJ,|l3) : Semi- 
classically, vertex corrections describe the interference of 
two (time-reversed) electrons undergoing an interaction 
process (see Fig.EJa)). Due to the transferral of a certain 
energy AE, the two electrons subsequently accumulate a 
phase difference, e lAEt , where t ~ t v . Therefore, ver- 
tex corrections with AE > 1/t v (characteristic for our 



3 



problem) do not effectively contribute to WL. (In con- 
trast, Coulomb interactions in d < 2 are dominated by 
low-energy transfers and vertex corrections are required 
to regularize infrared singularities |ll|.l 

We next turn to the discussion of the second family of 
corrections potentially altering our above results, corre- 
lated disorder-interaction processes as shown in Fig.^c). 
A preliminary indication as to the relevance of such con- 
tributions may be obtained by estimating the sample-to- 
sample fluctuations [18j of the Kondo temperature, Tk, 



7/2 , 
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(7) 



Substituting the results for the fluctuations of the local 
density of states 5v in weakly disordered d-dimensional 
metals Q, we obtain 



STk 



i i 

i i 

k F l (Jk)3 
1 1 

(fc F ;) 2 TJW 



tTk 



in (quasi) d = 1, 

in d = 2, (8) 

in d = 3, 



where L\ is the cross-section of a quasi one-dimensional 
wire. We assume, henceforth, that k-pl is sufficiently 
large, so that 5Tk -C Tk. While this condition seems 
restrictive in quasi-ld, it turns out to be always met in 
the WL regime, Atr WL <C (j Drude , realized in experiments 
(and assumed in this paper). 

More formally the role of correlations disor- 
der/interactions may be explored in terms of the dia- 
grams shown in Fig. Ht c )- On the face of it, these di- 
agrams are smaller by factors of l/(kpl) than the lead- 
ing contributions considered above (as quantum interfer- 
ence maintained across the impurity limits the momen- 
tum exchanged to values < l^ 1 much smaller than /cf). 
However, for very low T the enhanced infrared singu- 
larity caused by the presence of extra diffusion modes 
may over-compensate this phase space suppression fac- 
tor. Using that for T <C Tk, the bare interaction may be 
described by Fermi liquid theory , we find that only in 
one dimension the additional diagrams (oc T^ d+2 ^ 2 + T 2 ) 
lead to contributions of anomalously strong singularity. 
Specifically, we obtain a correction to the dephasing rate, 
1/%, C ~ n s T 3 / 2 /[vT£^(k F L ± ) 2 }. Therefore, for 



T < 



(fcpLj 



-Tk < T K 



(9) 



the separation disorder/interactions used above becomes 
invalid in d = 1. 

At very low T, yet another type of correction begins 
to play a role: The diluted Kondo impurities become in- 
distinguishable from a conventional disordered Fermi liq- 
uid with short-range momentum-conserving interactions 
|20| and the dephasin g ra te is determined by Altshuler- 
Aronov-Khmelnitsky [ill l2(i| type processes which, in 



our context, are encapsulated in the third family of dia- 
grams shown in Fig. Old) ■ Contributing only at order n| 
they generate corrections scaling as T 2 / 3 , T, and T 3 / 2 in 
d = 1,2, 3, respectively. Evaluating the prefactors (again 
in Fermi liquid theory) we find that these contributions 
become sizeable at temperatures 



(fc F 
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tT 2 
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-T K 

1 f ns 

tT k ) 1 /4 y„T K 



1/4 



rf = 3, 
d = 2, 

T K d=l. 



(10) 



In all cases the crossover scale is well below Tk and may 
arguably be neglected in all relevant experiments. Fur- 
ther corrections to 1 /t v of order u| and higher arise from 
clusters of two and more magnetic impurities which are 
sufficiently close that the inter-impurity coupling domi- 
nates over the Kondo effect jilj . 

Fitting to Aa = 2r 2 a Dr udo/ d d q(Dq 2 + ^(T))- 1 , 
we next relate our results to the T-dependent dephasing 
rates, r ( ^ 1 (T) extracted from experiments. Comparison 
with the energy-resolved representation J3J results in 



Jdef'(e)r v (e,T)^ 
exp[J def(e) In r v (e,T) 
-Jdef(e)/r v (e,T) 



d = 1,3, 
d = 2, 
lobt v > 1, 



where the last line applies to the case where a strong 
magnetic field is present and the Cooperon can be ex- 
panded in 1/ (u>bt v ). (ojb is the "cyclotron" frequency of 
the Cooperon.) 

Just two parameters enter the expression for 1/r^: the 
Kondo temperature Tk and the dimensionless density 
of magnetic impurities tis/(^Tk) <C 1. The T-matrix 
in Eq. 0] can be evaluated for arbitrary values of T/TW 
and e/Tk using the numerical renormalization group |22( 
generalized to the calculation of dynamical quantities in 
Refs. I23I Eif ). The result is shown in Fig. [21 where we 
define 25| and determine the Kondo temperature Tk = 



t£ 0) from the T = susceptibility \ = (.g^s) 2 /( 4T K 
For T 3> Tk, we obtain (correcting factors 2 in Q and 
16 in 14]) 1/t v « 37ms /(8u In 2 [T/Tk] ) giving rise to a 
very weak T dependence for T > T K - For T < 0.3T K 
one observes a crossover regime where 1/t v varies almost 
linear in T while at lowest T one obtains the expected 
T 2 behavior. This low-temperature regime can also be 
calculated analytically using Fermi liquid theory [l9| . 
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with c FL ~ 0.927,0.946,0.969 in d = 1,2,3 and c FL = 1 
for u!bt v ^S> 1 (correcting prefactors in Ref. 0). Acci- 
dentally, the dependence of both the analytical and nu- 
merical results on ujbt v (and on d) is very weak, implying 
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FIG. 3: Universal dephasing rate calculated via NRG. Acci- 
dentally, there is almost no dependence on the dimension d. 
The dot-dashed line shows the result of a Nagaoka-Suhl re- 
summation 26) , — = t? 2 — •>„ ,7, , using our definition 
of Tk 25]. While the Fermi liquid behavior is recovered for 
T < O.lTk, the calculated prefactor of the T 2 behavior (solid 
line) is not exactly reproduced due to numerical problems (23| . 



that the usual fits of the magnetoresistivity can be used 
to determine 1/t v . 

In the comparison to concrete experiments one needs 
to account for the interplay of dephasing due to magnetic 
impurities and due to Coulomb interactions ^]J. Since 
the latter are controlled by infrared divergences in d < 2, 
the respective rates do not add. Instead one needs to 
solve, e.g. in quasi one-dimensional systems, the equation 



1 




2/(3v,s) 



(12) 



where the first term describes the self-consistently calcu- 
lated effects of Coulomb interactions while l/r Vi s is the 
dephasing rate due to the magnetic impurities. The first 
(second) line holds when the Coulomb dephasing (Kondo 
dephasing) dominates. 

In conclusion, we have shown that the dephasing rate 
due to diluted magnetic impurities can be calculated di- 
rectly from the inelastic cross section (0J introduced in 
fl5j| . This result is valid for all types of diluted dynami- 
cal impurities as long as typical energy transfers remain 
larger than 1/t v . In the case of Kondo impurities the de- 
phasing rate depends on just two parameters, the Kondo 
temperature Tk and the dimensionless density of mag- 
netic impurities ns/(^Xk) € l A measurement of t v 
for spin-1/2 impurities along with an independent exper- 
imental determination of Tk and ng would put our theory 
to a parameter-free test [up to fitting the Coulomb back- 
ground, Eq. (|12f) ]. An open question is how spin-orbit 
interactions and disorder influence the Kondo effect in 
systems with larger spins, e.g. Fe in Au 0. 
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